Ground State of Electrons in Bilayer Graphene 
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Based on the four-band continuum model, we study the petaloid-loop-current (PLC) state for 
electrons in bilayer graphene at the charge neutrality point. The present work resolves the puzzles 
that (a) the energy gap increases significantly with increasing the magnetic field B, (b) the energy 
gap can be closed by the external electric field of either polarization, and (c) the particle-hole 
spectrum is asymmetric in the presence of B, all these as observed by the experiment. We also 
present the prediction of the hysteresis energy gap behavior with varying B, which explains the 
existing experimental observation on the electric conductance at weak B. The large energy gap of 
the PLC state is shown to originate from the disappearance of Landau levels of n = and 1 states 
in conduction band. By comparing with the existing models and the experiments, we conclude that 
the PLC state is a possible ground state of electrons in bilayer graphene. 

PACS numbers: 73.22.Pr,71.70.Di,71.10.-w,71.27.+a 
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The study of bilayer graphene (BLG) is a focused area 
in the condensed-matter physics because of the poten- 
tial application of BLG to new electronic devices [J- 
4]. One of the fundamental subjects is to explore the 
physics of the ground state of electrons in BLG. A num- 
ber of experiments performed on high quality sus- 
pended BLG samples have provided the evidence that 
the ground state is gapped at the charge neutrality 
point. In particular, a recent experiment by Velasco et 
al. 0] has observed that (i) the ground state is insu- 
lating in the absence of external electric and magnetic 
fields, with a gap E gap ps 2 meV that can be closed 
by a perpendicular electric field of cither polarization; 
(ii) the gap grows with increasing magnetic field B as 
-Egap = A + V 'a 2 B 2 + Aq with A s» 1 meV and a ps 
5.5 meVT -1 ; and (iii) the state is particle-hole asymmet- 
ric. On the other hand, theories have predicted various 
gapped states, such as a ferroelectric-layer asymmetric 
state [§-[l3[ or quantum valley Hall state (QVH) a 
layer-polarized antiferromagnetic state (AF) 15[ ; a quan- 
tum anomalous Hall state (QAH) [l2l . [la . fl7j j , a quantum 
spin Hall state (QSH) [l2|, LL3]> an d a superconducting 
state in coexistence with antiferromagnetism (SAF) fl8j j . 
The ferroelectric-layer asymmetric and QAH and QSH 
states all have been ruled out by the experiment Q . The 
SAF state is excluded because the real system is insu- 
lator. The AF state cannot reproduce the gap behavior 
with varying the magnetic field. A candidate proposed 
very recently is the loop-current state that is a broken in- 
version and time-reversal symmetry state 19|, 2^]. This 
state has been studied by numerical diagonalization of an 
effective mean-field Hamiltonian for a finite size lattice 
in presence of the magnetic field [l^] and by a naly tical 
solving a two-band continuum model (2BCM) [20(. By 
viewing the structure of the current loops connected to 
each lattice site [2(j, we here name it as the petaloid- 
loop current (PLC) state. Whether the model of PLC 
state agrees with the experimental observations on the 



electronic properties of BLG remains a question. 

In this work, using the four-band continuum model 
(4BCM) for electrons with finite-range repulsive inter- 
actions in BLG, we study the PLC state at the charge 
neutrality point and compare the results with the exper- 
imental observations. We investigate the gap behavior 
of the PLC state with varying the magnetic field B, and 
the particle-hole asymmetry spectra at finite B, and the 
phase transitions in the electron system in the presence 
of the electric and magnetic fields. We will show that the 
puzzles (i)-(iii) of the experimental observations can be 
resolved by the present model. 

The 4BCM. The top view of the lattice structure of a 
BLG is shown in Fig. 1(a). The unit cell of BLG contains 
four atoms denoted as ai and bi on top layer, and a2 
and b2 on bottom layer. The lattice constant defined as 
the distance between the nearest-neighbor (NN) atoms 
of a sublattice is a ps 2.4 A . The distance between the 
two layers is d ps 3.34 A . The energy of intralayer NN 
[between ai (&2) and bi (b^)] and interlayer NN (between 
bi and a.2) electron hopping are t ~ 2.8 eV and t\ ps 0.39 
eV, respectively. The reciprocal lattice of a single layer 
and the first Brillouin zone are shown in Fig. 1(b). 

For the carrier concentration close to the charge neu- 
trality point, we need to consider only the states close 
to the zero energy. Since the momenta of the low-energy 
states are close to the Dirac points K = (47r/3, 0) and 
K' = —K, we here define the operator 



pf _ ( J t t t \ 



where v — K or K 1 



creates a spin-tr electron of 



momentum k in valley v of I sublattice, and k is mea- 
sured from the Dirac point K (K 1 ) and confined to a 
circle k < 1/a of K (K 1 ) valley as shown in Fig. 1(b). 
In the momentum space, the Hamiltonian describing the 
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FIG. 1: (Color online) (a) Top view of the bilayer graphene. 
Atoms ai (a2) and bi (b2) are on the top (bottom) layer, (b) 
The first Brillouin zone (hexagon) and the two valleys K and 
K' in the reciprocal-lattice space of honeycomb lattice. 
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FIG. 2: (Color online) Some of the current loops connected to 
the black site on a sublattice. The sign factor ± in Iu{j — i) = 
±Ji(\j — i\) is + (-) for the electron motion from the black 
site i to the red (white) site j. 



noninteracting electrons is given by 
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where e„& — eo(s v k x + ik v ), s v — 1 (-1) for k in the valley 
K (K') : and e = V3t/2. We hereafter use the units of 
eo = l and a = 1. 

The interaction Hamiltonian is given as 

H' = U ^ SnttfSriiji + - ^ vu^j5n u 8n Vj (3) 



U 



where 8nu a is the number deviation of electrons with 
spin a from its average occupation at site i of sublattice 
I (hereafter denoted as li for short), Snu = Snu-f + Snui, 
U is the on-site interaction, and vu.i'j is the interaction 
between electrons at sites li and I'j. Within the mean- 
field approximation (MFA), since the interaction vnj>j 
appears in the exchange self-energy, it can be considered 
as a finite-range interaction by taking into account the 
screening effect due to the electronic charge fluctuations 
2l|. The total Hamiltonian Hq + H 1 satisfies the particle- 
hole symmetry |2l| . 

The PLC state. Under the MFA, the self-energy £«/ (k) 
of electrons in the PLC state is obtained as 



h'{k) = ^] vu,i'j (QtgCj,^; exp 



) exp(zfc • d) (4) 



with (c u<7 c\ lja 



Rw{d) + ilui {a), where d is a vector 

from the position li to I'j. The quantity (d) is propor- 
tional to a current. First, we consider the diagonal self- 
energy £j;(fc) for the same sublattice I. The formulae for 



Ru(d) and Iu(d) are derived in the supplement. As men- 
tioned above, only the low-energy states with momenta 
close to the Dirac points need to be considered. Then by 
expanding the self-energy with respect to k in the two 
valleys and taking only the leading terms S;/(±if), we 
obtain (see the supplement) 

T,u(±K) = nv c ±div s 



n 



di 



k 

lj H (( C \Ki+ka C lK'+k<y) - (4K+ka C lK+ka)) (5) 



N 

k 

v c = v(d) cos 2 (i? ■ d) 
v s = ^u(d)sin 2 (i? • d) 

d^O 

where n and di come respectively from the parts Ru(d) 
and Iu(d), the k summations here are confined to a single 
valley, v(d) — Vnjj, and is the total number of unit 
cells on single layer graphene. The quantity ri can be 
written as r; = — Si/ 2 with 61 as the average electron dop- 
ing concentration on sublattice I. For the homogeneous 
system at the charge neutrality point, we have r; = 0. 
The quantity div s is the order parameter of the PLC 
state. In Fig. 2, we draw out all of the current loops on 
the same sublattice connected to a given site i. Clearly, 
the total current density at site i is zero. The existence 
of the current loops breaks the time-reversal symmetry. 

Second, when the quantity (a) with I ^ V is not 
vanishing, it breaks the translational invariance of the 
system. This case happens only when the interactions are 
strong enough. For weak to medium interactions, we here 
assume all the currents between the bonds of different 
sublattices are negligible small (see the supplement and 
Ref. [ljj. 
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A particular property of the PLC state is that the 
order parameter term is odd for interchanging the val- 
ley index. We suppose d ai v s — —d}, 2 v s = — A x and 
db t v s — —d a2 v s = A 2 that means the breaking of the 
layer inversion symmetry. As a result, the effective 
MFA Hamiltonian H v k is obtained by adding the di- 
agonal matrix Diag(— s„Ai, s v A 2 , —s v A 2 , s„Ai) to H° k . 
Note that the matrices H v k and H_ v k are related by 
H v k — SH- v ^kS, where S is a 4 x 4 matrix with ele- 



ments S 



fJLV 



If tpu is an eigenfunction of H, 



with eigenvalue E£ (/x = 1,2,3,4), then Sip% is an eigen- 
function of H- v -k with the same eigenvalue. Therefore, 
only the eigenstates in a single valley need to be solved, 
which give rise to the whole energy spectra. 

1. The PLC state at B = 0. With the wave functions, 
the order parameters arc determined by 



Ai = 



A 2 = 



V3v s 
2V 

2V 



k (J 
k f! 



I^T) 



(6) 
(7) 



where / is the Fermi distribution function, ipj* is the z/th 
component of the eigenfunction ip?, V = \/3N/2 is the 
total area of the one layer lattice, and the k summation 
runs over a single valley. From the 2BCM, we know that 
the valence and conduction bands are connected to the 
electronic motions in the ai and b 2 sublattices. There- 
fore, the energy gap between the valence and conduction 
bands is determined by 2Ai. To reproduce the experi- 
mental data | Ai| = 1 meV at the charge neutrality point, 
v a needs to be 5.8eo. Supposing the effective interaction 
v(r) Rj e 2 exp(— qor)/re with e rj 3 as the screening con- 
stant of high frequency limit of BLG, we can obtain the 
desired value v s with qo sa 4.54/a. Another coupling 
constant is obtained as v c 1.94eo. 

2. The PLC state at finite B. In the presence of the 
magnetic field B applied perpendicularly to the sample 
plane, we take the Landau gauge for the vector potential, 
A = (0, Bx). With this gauge, the y component momen- 
tum k y is a good quantum number. Replacing the vari- 
able x and the operator k x = — i\J x with the raising and 
lowering operators <r and a, k y +Bx — y5/2(a'+fl) and 



i^Bj2(c 



t. 



a) , we can rewrite the effective Hamil- 



tonian in real space. The eigenfunction ipKn(j) ( w hh A 1 
= 1,2,3,4) can be expressed as 



IX 



2 fj 



Kn 



V 



X K., ■ 
X Kn ( t>n- 
-ix Kn 4> n 



-2(j)J 



(8) 



for n > 2, where (f> n (J) * s the n ^h level wave func- 
tion of a harmonic oscillator of frequency 2B and mass 
1/2 centered at x c = —k y /B. The vector X^ Cn — 
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FIG. 3: (Color online) (a) Landau levels Ej Cn in the valence 
(squares) and conduction (circles) bands at B = 1 T. The 
solid and dashed lines represent continuum conduction and 
valence energy bands, respectively, at B — with the abscissa 
being momentum k. (b) The energy gap i?gap (diamonds) as 
function of the magnetic field B compared with the AF result 
and the experimental observation (Exp, Ref. 0). 



(^Kri'^iYn'^/fri'^ifn)* (with the superscript t meaning 
the transpose of the vector) and the eigenenergy fit 
are determined by 



with 



H 



K r. 
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V2Bn A 2 
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The vector is normalized to unity. For each n > 2, 
the four energy levels appear at the valence, conduction, 
and other two bands about ±ii far from the zero energy, 
respectively. For n = 1, there are only three states with 
and the other three components and eigenvalues 



An _ 

V K1 



are determined by the upper left 3x3 matrix of Hxi- 
For n — 0, we have only one state X^ = (1, 0, 0, 0) and 
EjfQ = — Ai. To self-consistently determine the order 
parameters, we need to calculate the quantity d;. The k 
summation in Eq. ([3]) now is changed to the sum over 
the quasiparticle states. The order parameters are now 
determined by 

Al = ^E/(^n)[l^nl a -|^J 9 ]. (10) 



4tt 



a ft 



A 2 = 



VZvsB 
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(ii) 



Here B/2tt is the degeneracy of a Landau level, and B is 
in the unit of B = hc/ea 2 = 1.105 x 10 4 T. 

The solution to the Landau levels at B = 1 T is shown 
in Fig. 3(a). Only the levels in the conduction and va- 
lence bands are depicted. For n — there is only one 
state of energy —Ai(B) in the valence band. For n = 1, 
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FIG. 4: (Color online) The hysteresis curves for Ai(B), 
A 2 (B) (left) and £ gap (B)Ai(0) (right). 



there is a level E v Kl slightly above —Ai(B) in the valence 
band. There is no state in the conduction band for n = 
and 1. Only when n > 2, the level E c Kn in the conduction 
band appears. The energy gap is 



-Egap — E C K2 — E Kl . 



(12) 



Clearly, the particle-hole symmetry is no longer valid at 
finite B, in agreement with the experiment [7|. Fig. 3(b) 
shows the energy gap E gap of the PLC state as function 
of the magnetic field B. The AF calculation of the same 
4BCM and experimental results for E gap are also plotted 
for comparison. Here, the only fitting parameter is v s 
for reproducing Ai(0) = A at B = 0. The theoretical 
result for E gap of the PLC state as a function of B is in 
surprisingly good agreement with the experiment 0- 

The above solution to the order parameters is only in 
the branch of A x > and A 2 > 0. At weak magnetic 
field B > 0, there is another branch in which Ai < and 
A2 < 0. In this case, the two levels of n = and 1 appear 
in the conduction band but not in valence the band, and 
the energy band is given by E gap = Ej C1 — E V K2 . In 
Fig. 4, we show the hysteresis curves for the PLC order 
parameters Ai(B) and A 2 (B) and the energy gap E gap in 
the range -0.5 T < B < 0.5 T. For \B\ < 0.18 T, there are 
two branches for E gap . In the lower gap branch, the gap 
decreases with increasing \B\. This behavior of E gap is 
in qualitative agreement with the experiment observation 
by Weitz et al. [5| that indicates two peaks in the electric 
conductance appearing at B p ±0.04 T (where the real 
gap reaches the minimum), respectively. 

The AF state. In the AF state, the two layers have 
opposite magnetizations that stem from the on-site in- 
teraction U. We denote the magnetizations in a unit cell 
as (m ai , m,b 1 , m a2 , mb 2 ). The order parameters are given 
by -U(m ai ,m bl ,m a2 ,m b2 ) = (-Ai, A 2 , -A 2 , Ai). 

At B = 0, the equations for determining the order 
parameters happen to be the same as Eqs. (|6]) and (JT]). 
To reproducing the energy gap of the experimental data, 
U = v s sa 5.8eo is needed. (See supplement.) 

At finite B, the spectra of the electrons are different 
from the PLC state. At valley K, the Landau levels have 



the similar distribution as the PLC state. In particular, 
there are no levels of n — and 1 in the conduction 
(valence) band if Ai > (<)0. However, at valley K', 
these levels appear but their counterparts disappear in 
the valence (conduction) band. As a result, the energy 
gap is not as large as the PLC state but is 2Ai(B). The 
self-consistent solution for the energy gap in the presence 
of B is shown in Fig. 3(b). Obviously, the B-dependent 
AF gap cannot explain the experimental result. 

The PLC state under external electric field. When an 
external electric field E is applied perpendicularly to the 
BLG plane, there is a potential difference 2u — Eed/e 
between the two layers. The Hamiltonian H v k now is 
obtained by adding the diagonal matrix Diag(u + r\v c — 
s v A 1 ,u + r 2 v c + s v A 2l -u-r 2 v c -s v A 2 , -u-nVc + SyAi) 
to H® k . Here the terms r;t> c appear because of the electric 
polarization by E. Note that n has the same sign of u, 
and thereby H v k{u) = SH_ v -k{—u)S, which means that 
the order parameters are even functions of u. The model 
shows that if an external electric field in positive direction 
closes the energy gap, then an electric field in negative 
direction does it either. For the sake of illustration, we 
here consider the case of B > and u > 0. The results 
for other cases can be deduced by the symmetry of the 
Hamiltonian. For B > and u > 0, we still have two 
cases: Ai > and Ai < 0. Here, we consider the case 
of Ai > 0. The discussion can be extended to the case 
of Ai < 0. At B = 0, the effective gap is u + r\v c — 
A\. The critical potential uq closing the effective gap is 
obtained as uq « 0.48Ao ~ 0.48 meV. The critical field 
of the experimental data 0] is E w 1.25 mVA -1 , which 
corresponds to uq sa 0.69 meV (using e 3). The theory 
is in qualitative agreement with the experiment. 

Since the system satisfies the particle-hole symmetry 
at B = 0, we can take the chemical potential as zero 
for the system at the charge neutrality point. Then, the 
level u + r\v c — A\ is occupied if its sign is negative, 
otherwise it is empty. Therefore, with increasing u from 
0, the system undergoes a phase transition at u = uq 
from the state with the level u + r\v c — A\ occupied 
to the state with the level empty. Thus, to search the 
critical u where the gap closes at finite B, we study the 
phase transition. 

At finite B, a state at (B, u) can be obtained by contin- 
uously changing the parameters B and u from the state 
at (0, Ui). If Ui > uo(0), then the level u + r\v c — A\ is 
empty. Note that u + r\v c — A\ is the only Landau level 
Eku of n — at finite B and there is another level of 
n = 1 close to it similarly as the case of u = 0. So the 
two levels of n = and 1 in the K valley keep empty 
on the path from (0,Uj) to (B,u). On the other hand, 
if one starts from an initial state with Ui < uq(0), then 
the two levels of n = and 1 keep filled. (We denote 
the filling number as Jk — and 1 for the two levels 
empty and filled, respectively.) We thus have two states 
at (B,u). To judge which one is the real ground state, 
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support Q for the prediction for the hysteresis energy 
gap behavior with varying B. These facts show that the 
PLC state is a possible ground state of electrons in BLG. 
The model explored here can be useful for understanding 
the physics of the electrons in BLG that is expected as a 
new generation of semiconductor. 

This work was supported by the National Basic 
Research 973 Program of China under Grants No. 
2011CB932702 and No. 2012CB932302, NSFC under 
Grant No. 10834011, and the Robert A. Welch Foun- 
dation under Grant No. E-1146. 



FIG. 5: (Color online) Phase boundary uo(B) between the 
two phases (fx, fx 1 ) = (0, 1) and (1, 1) (red solid line). The 
solid points and the diamonds are converted from the critical 
E of the experimental data (Ref. 0) using e ~ 2.24 and 3, 
respectively. The inset shows the result for Ai < in the 
range < B < 0.18 T. 



we compare their ground-state energies. With the energy 
judgment, the ground state at (B,u) is uniquely deter- 
mined. At the critical potential uq(B), the two states 
have the same ground-state energy. 

In Fig. 5, we exhibit the result for uq(B) as function 
of B and compare it with the experimental data Q • The 
experimental data are obtained by converting the criti- 
cal electric field E to uq according to uo — Eed/2e with 
the dielectric constant e « 2.24 (solid points) and 3 (di- 
amonds). As seen from Fig. 5, the behavior of uq(B) 
by the theoretical calculation is in fairly good agreement 
with the experiment [7J with e w 2.24. 

As already seen, there is another solution of A x < in 
the range < B < 0.18 T. We show in the insert in Fig. 
5 the phase boundary for this case. We see that the state 
of Ai < in 0.07 < B < 0.18 T is unstable with respect 
to a small E. The range for the stable state of Ai < 
is reduced to \B\ < 0.07 T, with [B max \ = 0.07 T close 
toward to the experimental data [5[ \B P \ = 0.04 T. 

Summary. With the MFA to the 4BCM, we have stud- 
ied the PLC state of the electrons with finite-range repul- 
sive interactions in BLG at the charge neutrality point. 
Because of the broken inversion and time-reversal sym- 
metry in the PLC state, the Landau levels of n = and 1 
states in the conduction/ valence band disappear. With 
only one coupling constant v s fitting the experimental 
gap Ao at B — 0, the obtained energy gap at finite B is 
in surprisingly good agreement with experimental data 
[3] • The results for the phase transition in the system in 
the presence of external electric and magnetic fields, and 
the particle-hole asymmetry spectra in the presence of B 
are also in good agreements with the experimental obser- 
vations Q- There is also the intermediate experimental 
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